DO-TH-99/23 
December 1999 



Renormalization of the nonequilibrium 
dynamics of fermions in a flat FRW universe. 

Jiirgen Baacke f] and Carsten Patzold Q 



Institut fiir Physik, Universitat Dortmund 
D - 44221 Dortmund , Germany 



Abstract 

We derive the renormalized equations of motion and the renormalized energy- 
momentum tensor for fermions coupled to a spatially homogeneous scalar 
field ( inflaton) in a flat FRW geometry The fermion back reaction to the 
metric and to the inflaton field is formulated in one-loop approximation. 
Having determined the infinite counter terms in an MS scheme we formulate 
the finite terms in a form suitable for numerical computation. We comment 
on the trace anomaly which is inferred from the standard analysis. We also 
address the problem of initial singularities and determine the Bogoliubov 
transformation by which they are removed. 
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1 Introduction 



The production of fermions during inflation |]J was for a long time consid- 
ered to be less interesting than the one of scalar quanta. The production of 
scalars coupled to the inflaton field is characterized, during preheating [FJ, |], 
by parametric resonance bands 0, [|], a feature that leads to an abundant 
production of such quanta. Numerous authors have investigated various as- 
pects of the production of scalar quanta, including the back reaction to the 
inflaton field and to the scale parameter [5 — 14]. Also, parametric reso- 
nance plays an important role in the analysis of the late time behavior of 
these coupled systems |TE| , |T7|] , in the large- N limit. For fermions a true 



parametric resonance cannot develop, as the fermion number is limited by 
the Pauli principle. So fermions were not considered in the inflationary or 
preheating stage of inflation but thought to be produced afterwards via the 
decay of scalar fields. Recently it was observed |18|, |19|] , that the production 
of fermions during the preheating stage, via their coupling to the inflaton, is 
characterized by resonance-type bands within which the occupation number 
comes near saturation. On the other hand, the production of fermions via 
decay of a boson into an fermion- ant if ermion pair leads to a sharp momen- 
tum spectrum. As a consequence the coherent production of fermions in the 
inflaton background may largely exceed the production rate obtained via the 
decay of scalar fields. 

This phenomenon has lead several groups to (re) consider the production 
of fermionic degrees of freedom during the inflationary stage. The nonper- 
turbative production of fermions in an axial background field was consid- 
ered in |2(J . Very massive fermions may produce characteristic spikes in the 
spectrum of primordial perturbations, leading to observable features in the 



measured CMB anisotropy and in the power spectrum [21]. In supergravity 



gravitinos will necessarily be a part of the theory, their parameters can be 
related to the scale of super symmetry breaking. The coherent, non thermal 
production of gravitinos has been considered in Refs. [22, |23], . It may lead 



to an exceedingly low reheating temperature and then becomes a problem for 
many models of the early universe [^, p6j] . This issue still is still open and 
being considered by several authors. One of the problems is the modification 
of the parametric resonance bands by the evolution of the scale parameter, 
or by an evolution of the inflaton field modified by the back reaction. 

Here we will consider the renormalization of the fermionic back-reaction 
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to the inflaton field and to the energy momentum tensor in an FRW geom- 
etry, in one- loop approximation. Though this seems a more formal aspect, 
it is important for a reliable numerical computation of these back-reaction 
effects. As long as one just considers the number of produced particles this 
is not necessary as the adiabatic particle number is finite per se. For the 
computation of the other quantities simple normal ordering is not sufficient 
and leaves at least logarithmic divergences. 

The renormalization of nonequilibrium dynamics of fermions in Minkowski 



space has been considered by us previously JjJ|, but this analysis has to be 
adapted to the FRW geometry. In particular dimensional regularization is 
modified by the time dependent scale factor. While the flat space counter 
terms still renormalize the divergent loops the finite terms are different. Fur- 
thermore, the structure of the energy momentum tensor in the FRW geom- 
etry is modified as compared to the one in Minkowski space. So several new 
renormalization constants and finite parts have to be determined. For the 
case of scalar quantum fluctuations we have already presented [§, [K| our 



scheme JTTJ for renormalizing the equation of motion for the inflaton and 
the Friedmann equations. The methods used there, well suited for numerical 
computations, can be adapted to the case fermions. In this sense the present 
publication merges our previous work on fermion systems in Minkowski, and 
scalar fields in FRW space. Another viable scheme is the well-known adi- 



abatic subtraction ||27|| . While this has been worked out and numerically 
implemented for the cosmological nonequilibrium evolution of scalar fields 
in Ref. f5|, W^, an application to fermions has not yet been adapted for the 
requirements of numerical computation. 

Of course not all the results presented in this publication are new. The 
renormalization of the energy momentum tensor of free fermions in a FRW 



background has been considered long ago, see e.g. [£8], |29|, [!0|; other regu- 
larization schemes have been discussed, e.g., in Refs. |[T, 32]; the one-loop 
renormalization of a Yukawa theory is of course standard. However, while 
the structure of the renormalization constants is known, we have to rederive 
them here within the nonequilibrium formalism. An important point is their 
independence of the initial conditions. Furthermore, the remaining finite 
parts have to be written in a form suitable for numerical computation. As 
usual, the formulation of the finite parts is the most cumbersome part of 
renormalization, and this has not been discussed before. 

The plan of the paper is as follows: in section 2 we recall the basic 
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formulae for a flat Friedmann universe, with 3+e space dimensions. In section 
3 we embed into this geometry the quantum field theory describing massive 
fermions, Yukawa-coupled to the scalar inflaton field. The energy momentum 
tensor for the classical field and the fermionic quantum fluctuations is given 
in section 4. In section 5 we determine the renormalization constants, and 
thereby the precise formulae for the finite parts. We also include the trace 
anomaly that cannot be derived from considering the FRW metric alone. We 
briefly digress on the question of initial conditions in section 6, conclusions 
are presented in section 7. 



2 FRW cosmology 

We consider the Friedmann- Robertson- Walker metric with curvature param- 
eter k — 0, i. e. a spatially isotropic and flat space-time. The line-element 
is given in this case by 

ds 2 = dt 2 -a 2 (t)dx 2 . (2.1) 

The time evolution of the cosmic scale factor a(t) is governed by Einstein's 
field equation 

(I + 5Z)G txv + 5a ^ + 5(3 (2) # llv + 8 1 H ltv + 5Ag ^ = -k{T^) (2.2) 

with k = 8nG. 

The Einstein curvature tensor G^ v is given by 

G^ = R^v — \g^ v R ■ (2-3) 

The Ricci tensor and the Ricci scalar are defined as 

R tw = R lv\ ' ( 2 - 4 ) 

R = gTRp,, (2.5) 

where 



3 



In FRW geometry with n — 1 space dimensions, as needed for dimensional 
regularization, the relevant quantities are given by 



Ru = (n-1)- 
a 



/ \ 2 

i2{J = i2 = 2(n-l)- + (n-l)(n-2) (- 

a \a, 



(2.7) 



It is customary to introduce the Hubble parameter a/ a and to use R = (n — 
l)(2H + nH 2 ). The tensors ^H^ u , ^'H^, and if^ arise from the variation 
of terms proportional to i? 2 , R al3 R a/ 3, and R 01 ^ 5 in the Hilbert-Einstein 
action. Their general definition is given, e.g., in pOfl . Though the fermion 
loop integrals will not produce divergences proportional to these tensors we 
will need them as they produce the trace anomaly. The explicit expressions 
in FRW geometry with n — 1 space dimensions (see, e.g. flOU ) are presented 
in Appendix A. 

The Friedmann equations are obtained by inserting the explicit expres- 
sions for the geometrical tensors and for the energy-momentum tensor into 
the Einstein equations. As this is done after renormalization we can set 
n = 4. Then 

3H 2 = n£ ren (2.8) 
R = K% en . (2.9) 

Here £ = T u and T = Tfi. For the numerical evolution it is sufficient 
to consider just the first equation. Then the second one is fulfilled as a 
consequence of covariant energy conservation 

E + H{AE-T) = . (2.10) 



3 Lagrangian and equations of motion 

We consider the production of fermions by a scalar field $ in a Heisenberg 
state in which the scalar field has a time-dependent, but spatially homoge- 
neous expectation value 

m = (<&(*)> (3.i) 



4 



The Lagrangian density of the model in Minkowski space is given by 

c = l±^d^d^ - v($) - ^-^m 2 - (c + sc)R$ - 



+ip (i^df, -m f -g®)ip (3.2) 



where 



+!l±^4 3 + (3.3) 

In addition to the scalar mass M, the quartic self-coupling A and the confor- 
mal coupling £ we have introduced a tadpole coupling a , a trilinear coupling 
k , and an additional coupling to the curvature tensor (RQ. For rrif ^ such 
terms are generated by the fermion loop and at least the counterterms are 
needed for the purpose of renormalization. Alternatively one may generate 
the fermion mass by a shift in $, starting with a massless fermion; then the 
divergent terms with odd powers of $ are generated from the terms with 
even powers of <£> obtained in the massless case. We have written the counter 
terms explicitly from the outset though their form is obvious. For the wave 
function renormalization this is necessary: if we use conformal time and con- 
formal rescaling of the fields in n = 4 — e dimensions, then part of the kinetic 
term, including the divergent wave function renormalization, reappears as 
part of the conformal coupling (see below). 

The action in curved space is obtained in the usual way by introducing 
covariant derivatives and the covariant measure, we have anticipated already 
the couplings to the curvature tensor. We will implement dimensional regu- 
larization in the way of using n — \ space components, so we will have to use 
FRW geometry with n — 1 = 3 — e space components, as well. 

The equation of motion of the condensate is given by 



[1 + SZ) 



+ (n- l)-4> 
a 



+ ^ + (C + K)R + (£ + 6t)R4>(t) (3.4) 



where the back reaction of the fermion field is given, in 1-loop and Hartree 
approximation by the term g(ipip), the back-reaction of the scalar field may 
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be included in one-loop, Hartree, or large-iV approximation. The renor- 
malization of the back-reaction of the scalar fluctuations has been discussed 



recently within our scheme in Refs. || |TD , and in the adiabatic subtraction 
scheme in Ref. |8], |15| . So we do not consider these fluctuations here. 
It is convenient to introduce conformal time 



dt = a{r)di 



and to rescale the field d> via 



r = a 



-l+e/2. 



r)0(r) 



The derivative with respect to conformal time will be denoted 

fir) = J-/(r) = d T f(r) . 
The classical equation of motion now takes the form 

(1 + 5Z)4>" + (C + 5()a 3 - £/2 R + a 2 

4 dV e ($) , X + 5X ~ 2 3- e /2/-T;\ r, 

where we have introduced the rescaled potential 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



H(0) 



(a + k)a- 1+£/2 ^ + 



M 2 + SM 2 



+ 



-2+e 72 



-4+2e 74 
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24 



(3.9) 



The covariant derivative in the Dirac equation is obtained using the n-bein 
formalism (see Appendix B). One obtains 



n — 1 a 



h°(d + ^~)il>+( fy fc V fc -m f - go?"* ) j, = . 



(3.10) 



The extra term in the first parenthesis can be eliminated by defining 
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furthermore, we introduce conformal time and obtain 

'id T -n(r)\ 4>(t,x) = . 



with 



H(t) = -iaV + \a(r)m f + #a e/2 (r)0(r)] (5 . 



(3.12) 



(3.13) 



Note that 7i is rescaled with respect to the Hamiltonian in Minkowski space 
with a factor a(r). So the Dirac equation takes its conventional form, with 
an effective time-dependent mass 



m(r) = a(r)m/ + ga e ^ 2 (r)4>{r) . 



(3.14) 



On account of the spatial homogeneity of the condensate field it is suitable 
to expand the Dirac field as 



(ftp 



(2tt) 3 2£ 



b PyS U PyS (r) + dl p ^ PyS (r) 



= +«p-x 



(3.15) 



with the time independent creation and annhiliation operators for quanta, 
whose mass mo = m(0) is determined by the initial state. Eq is the corre- 



sponding energy ym^ + p 2 . The creation and annihilation operators satisfy 
the standard anti-commutation relations 

{bp,s,bl >9 ,} = 2£o(27r) 3 5 3 (p-p')(W, (3.16) 
K, s ,4v) = 2£ (27r) 3 5 3 (p - p')8 ss/ . (3.17) 

For the positive and negative energy solutions we make the usual ansatz 



U p , s (t)=N id T + H p {r) f p (r) 



and 



Vp, s (r) = N [id T + n_ p (r)\g p (r) 
with the Fourier-transformed Hamiltonian 

H p (t) = ap + m{t)l3 . 



Xs 







Xs 



(3.18) 



(3.19) 



(3.20) 
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For the two-spinors Xs we use helicity eigenstates, i.e., 

pcrx ± = ± X ± ■ (3-21) 

The mode functions f p and g p depend only on p = |p|; they obey the second 
order differential equations 

d 2 



- im'(r) + p +m (r) 
d 2 

— + im'(r) +p 2 + m 2 (r) 
dr z 



f P (r) = 0, 
9p(t) = 0. 



(3.22) 
(3.23) 



If the fermion Fock space is based on the conformal vacuum state the modes 
start - in conformal time - as if the mass was independent of r for r < 0. 
The mode functions, which would be plane waves for t < 0, then satisfy the 
initial conditions 

/ p (0) = l , f p (0) = -iEo, (3.24) 
<fr(0) = l , g p (0)=tE , (3.25) 

so that g(r) = f*(r). The fermion condensate occurs in the equation of mo- 
tion as ga 3 ~ e / 2 (^tl)), the rescaling of the Dirac field changes this to ga e ^ 2 (ipip}. 
We define a fluctuation integral as 

d^p 



(U) = e / 



(2tt) 3 2£ 



V- p , s (t)V- p , s (t) 



, r d n -\ 

7 (27T) n - 1 



2En 



2p 2 



E + m 



\fp 



(3.26) 



(27r)- 1 2E 

We have written the mode integral in dimensionally regulated form. This 
corresponds to the usual prescription of introducing dimensional regulariza- 
tion 'after taking the Dirac traces'. Indeed the mode integral is a trace of a 
Green function. 



4 The energy-momentum tensor 

The energy-momentum tensor for a spatially isotropic background field and 
of the quantum fluctuations generated by this field is diagonal and of the 
form 

= di&g(£,V,V,V) . (4.1) 
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From the Lagrangian ( |3.2| ) we derive for the energy density of the condensate 
or background field 



£ C ond = ^^a- 4+e (0') 2 + V e (0) 



(4.2) 



+2(n-l] 



(n-2)(l + 5Z) 
An - 4 



-2+e 



Ha 



n-2 



-H 



2 12 



+2(C + 5()(n - l)a-^H4>> + b -^G tt + — , 

K K 

where V e (0) has been defined above, in Eq. ( |3.9| ). 

We have included the counter terms proportional to Gu and gu which ap- 
pear on the left hand side of the Einstein field equations. They are needed, 
as the other counter terms, in order to compensate the divergences arising 
in the fluctuation part of T a . The higher curvature counter terms on the left 
hand side of the Einstein equation ( |2.2|) are not needed for infinite renormal- 
izations, they will play a role, however, in the discussion of the trace anomaly 
at the end of section 5. 

Instead of the pressure we consider the trace of the energy-momentum 
tensor. We denote it by T = £ — {n — 1)V. For the background field we find 



cond 



+2(n-l) 



{n-2)(l + 5Z) 



4n - 4 



-2+e 



-1 J / 
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2 ~\ 2 
-H<p 



+2(n -!)(( + 5()a 



-3+e/2 



y, + (n-2)Ra> ~ 
n — 1 



6Z G „ U 5 A 
+ — ~ G u + n — 

K h K 



including again counter terms from the left hand side of the Einstein equa- 
tions. 

The fluctuation parts of the energy-momentum tensor are given, after 
conformal rescaling, by 



-4+e 



?/(2^ y -^K)^ 



(t) 



(4.3) 



2a 



-4+e 
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for the energy density and by 



r & (r) = ^— y(V>7P^> (4.4) 



a 



(27r) 3 2£, 

for the pressure. The fluctuation part of the energy-momentum tensor is 
given by 

(T£) fl = £-(n-l)P. (4.5) 

This results in 

%(r) = a -+.^|_J|_F_ p „(T) (0l r -1p)V- T „(r) 

=a ^!i^7^ {T)m(T)v -°'- (T) (4 - 6) 

So the trace can be computed by multiplying the fluctuation integral T by 
the time-dependent mass. 



5 Renormalization 

In order to develop the framework for renormalizing the one-loop equations, 
we write the equation of motion for the mode functions, Eq. (|3.22| ), in the 
form 

" ] f P (r) = -V(r)f p (r), (5.1) 



with 



V(t) = m 2 (r) - m(0) 2 - i [m'(r) - m'(0)] 



(5.2) 

Using the initial conditions ( |3.24| ) this equation can be recast into the form 
of an integral equation: 



f P (r) 



-iE(,r 



E Q Jo 



dr' S m[E (T-T , )]V(T , )f p (T'). 



(5.3) 
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Using this integral equation, the mode functions may be expanded with re- 
spect to the potential V(r). We split off the zeroth order (plane wave) 
contribution and an oscillating phase factor by writing 

f p (r) = e-^[l + h p (r)}. (5.4) 

The differential and integral equations satisfied by the function h p (r) can eas- 
ily be derived from Eqs. Q5.1J) and (|5.3|) , respectively. It may be decomposed 
as 

oo 

h P (r) = T,h p n) (r), (5.5) 



n=l 



where h p n \r) is of n'th order in V(r); we define further the inclusive sums 

oo 

h«) = £ h p m) . (5.6) 



m=n 



The truncated mode functions h p n \r) and h p n ^ satisfy a recursive set of dif- 
ferential and integral equations derived from Eqs. ( |5.1| ) and (|5.3[) ; these can 
be used to compute them numerically and analytically. Moreover the expan- 
sions (|5.5|) and ( |5.6|) are expansions with respect to 1/ E Q , we have used them 
[ |Tj|] in order to single out the leading divergent contributions in the various 
mode integrals. 

The integrand of the fluctuation integral J 7 , see Eq. Q3.26| ), can be written 

as 



i m o\ |f ^ |2 rn f m 



2Re h p (r) + \h p {r\ 



|2 



m(r) m"(r) m 3 (r) m{r)m^ 



E 4(£ ) 3 2(E )3 2(E )z 
m"(0) 

+ J7^i cos (2E t) + K F (p,T) . (5.7) 

The first terms on the right hand side lead to divergent or singular momen- 
tum integrals. The function K-p{r) can be considered being defined by this 
equation, i.e. by subtraction of the remainder of the right hand side from 
the left hand side. It behaves as (Eq)~ a and its momentum integral is finite. 
While other authors use the cutoff independence with scale dependent 



couplings in order to check their numerical scheme, we use this asymptotic 
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behavior as a numerical cross check. If defining K-p(r) by subtraction is 
numerically precarious, it can also be computed directly, using the reduced 

This is discussed in 



mode functions h p n \r) and 



For fermions these 



expressions get rather lengthy, we do not present them here. Using the trun- 
cated mode functions instead of subtracting plainly the divergent parts from 
the integrand constitutes another numerical and analytical cross-check. 
We decompose the fluctuation integral as 



F{t) = Fdivir) + ^singM + ^fin(T) 



(5.8) 



with 



div 



sing 



fin 



-2 
-2 
-2 



d n ~\> 

(27T)™- 1 



m(T) 



En 



m"(r) 



m 3 frl 



4(£ ) 3 2(E ) 



+ 



m(T)m 



o 



2(E ) 




4(£ 



, (5.9) 



m'(0) , s m"(0) , ,1 . 

7 sin(2£ r) + — ±L cos(2£ r) ,(5.10) 



(5.11) 



We have dropped dimensional regularization for the singular and finite parts. 
They are to be evaluated at n - 4 = e = 0. The divergent part has to be 
compensated by suitable counter terms. The part J sing displays a singularity 
in r at r = 0, a phenomenon related to the initial conditions. One gets 
rid of the problem by a Bogoliubov transformation |33|, 34 1, for fermions 
its explicit form has been derived in [[34], j35[. We derive the form of this 



Bogoliubov transform for the case under consideration in section 6. 

The fluctuation parts of the energy momentum tensor can be analyzed in 
a similar way. The integrand of the energy density 8$ can be expanded as 



m(r) 



-{E 



m J 



2Reh r 



\h v \ Im 

' P ' En 



m 2 (r) wig m /2 (r) m 4 (r) 
" °~^E~ + 2E + 8(Eor + 8{Eof 



mo m 2 (r)m§ 



8(E ) 



4(£ ) 3 



+ K E (p, t) . 



h'Jl 



h* 



m(r) 



(5.12) 
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Again Ke(p, t) is denned by this equation and it behaves as (Eq) 4 as E — > 
oo. There is no cosine term here and, therefore, no singular contribution. So 

£ fl (r) = £ fl ,div(T) + ffl,fl n (r) , (5.13) 



with 



Cdiv = 2a / — , : < — £j — h 771=- + N o + 



(5.14) 



(27T)"- 1 [ 2E 2E 8(£ ) 3 8(E ) 

m,Q m 2 (r)mg 
+ 8(Eof~ 4(E )3 j 

£*, = 2a- A J-0^K E (p,r). (5.15) 

The contribution of the fluctuations to the trace is particularly simple: as 
explained in the previous section, it is proportionla to the fluctuation integral. 
So, 7fl(r) can be decomposed as 

T fl (r) = T fl>div (r) + T fl>sing (r) + T fl;fin (r) , (5.16) 

where the three expressions on the right hand side are obtained from those 
of Eq. ( |5.8| ) by multiplying them with m(r)a~ 4+e . For discussing the renor- 
malization we need in particular the divergent part 

-4+ e f rf n_1 P Jm 2 (r) m(r)m"(r) m 4 (r) m 2 (r)mol 
Miv = « i^p 1 ! - ^ 4(£ ) 3 ~2(Eof + 2(E y J 

(5.17) 

The divergent terms J^div^div, and 7^ are proportional to local terms in 
4>(t) and its derivatives. These can be absorbed in the usual way by in- 
troducing the appropriate counter terms into the Lagrangian and into the 
energy-momentum tensor. 

The divergent parts of the fluctuation integral can be evaluated, e.g., 
using dimensional regularization. One finds 

m(r)ml 
with the abbreviation ^ 

1 f 2 4vr/i 2 a 2 (r) 1 
L ° = - + ln 2 - 7 • 5.19 



J- div = 2m"(r)L + Am 3 (r)L + , (5.18) 



3 This definition and the one of L and Lf below deviate from those in our previous 
work 1 18, 10| by inclusion of the factor 1/167T 2 . 
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The scale fi is introduced as usual by 'adjusting the dimension' of the various 
couplings, we have not displayed this rescaling explicitly. The initial mass 
only appears in the logarithm of L and in the last term of ( |5.18 ) which 



is finite. Therefore, as already found for the scalar fluctuations [[LI]], the 
dependence on the initial mass m = m(0) can be absorbed into finite terms, 
AZ, Act, .... In order to appreciate this independence of the initial mass it 
should be noticed that the individual divergent terms appear with prefactors 
up to mQ, and up to m\ in the energy momentum tensor. If one uses a three- 
dimensional 'computer' cutoff the necessary cancellations do not occur. 
We expand ( |5.18| ) using the explicit expression ( |3.14| ) for m(r). This leads 



to rather lengthy expressions, in particular we have to keep terms propor- 
tional to n — 4, as these become finite terms when multiplied with the 1/e 
of L . A further source of finite terms is the product L a e ^ 2 which appears 
when ga e l 2 T is evaluated in the equation of motion. We have used maple 
for computing the divergent terms and their finite remnants in the equation 
of motion. The second derivative term contains terms proportional to a", 
and therefore to the curvature scalar R. So the conformal couplings ( and 
£ are needed already in order to renormalize the equation of motion. In flat 
Minkowski space these couplings only appear in the pressure. 

Applying an MS prescription, the infinite renormalizations become 



sz 


= ~2g 2 L 


5a = 


-Agm 3 f L 


5M 2 


= -\2m)g 2 L 


5k, = 


-24m f g 3 


S\ 


= ~2^g 4 L 


s< = 


-\m f L 


5£ 


= ~\9 2 L 






L 


1 I 2 , 

16tt 2 \ e 


4-n-fi 2 
M 2 


-'■) 



(5.20) 



with 

1 f 9 Attn 2 

(5.21) 

We use a scale factor M instead of the fermion mass, as one may want to 
consider the case m./ = 0. The renormalization of the conformal coupling <5£ 
is due entirely to the appearance of 5Z in the conformal coupling term. 

The renormalized equation of motion is obtained by replacing the infinite 
parts of the counter terms by their finite remnants, and by including some 
additional terms. The finite remnants of the counter terms are 

AZ = 2g 2 L f Act = Agm)L f 
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AM 2 = l2m)g 2 L f An = 24m f g 3 L f 

AX = 24g 4 L f A( = -gm f L f (5.22) 



A£ = \g 2 L f 

with 

1 , M 2 a 2 (r) 
L f = —\n ±-L. 5.23 

These finite parts are time- dependent, the term \na(r) arises from the confor- 
mal scale factors via the product (2/e)a e//2 . Conceptually they are part of the 
remaining nonlocal finite fluctuation integrals, though they can be written 
in local form. 

The renormalized equation of motion, omitting the scalar back-reaction, 
reads 



(1 + AZ)f+ (C + A()a 3 R + a 
9 (2 



4 ^Hen(0) 



+ 16^ [^ a2mR + 4m o m + ^aHm - 2a 2 H 2 m^j (5.24) 
+9 (^sin g (T)+^ fi „(r))=0, 

where we have introduced the renormalized potential 

M 2 + AM 2 

V ren (0) = {a + Aa)a~ 1 <p+ a~ 2 2 

, K + AK o-o A + AA _ 4 74 /r orN 

+ g a 3 <f + 24 a 4 (5.25) 

that we will find again in the renormalized stress-energy tensor. We note that 
A£ and AZ have cancelled in the conformal coupling term; however, even 
if £ = 1/6 an additional conformal coupling oc g 2 is introduced by the first 
term in the parenthesis in Eq. ( |5.24| ). It does not vanish even for nif = 0. 
The cosmological constant counter term has not yet been fixed as it does not 
appear in the derivative of the potential. It will be determined below. While 
the equation of motion is now finite for r ^ it still contains ^ing(r) that 
is singular at r = 0. 
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The divergent parts of the energy density give, after dimensional regular- 
ization 

m 4 (0) m 2 (r)ml\ 



£ fl ,div = a" 4+e m'\r)L + m 4 (r)L + 



32vr 2 



+ 



8tt 2 



(5.26) 



We again expand this, using the explicit expression for m(r). In addition 
to the divergent terms found in the equation of motion we find terms that 
have to be compensated by the cosmological constant counter term and by 
the wave function renormalization of the gravitation field. This fixes these 
couplings to 



5k 



2 

mi 
3 ' 



(5.27) 
(5.28) 



again independent of the initial condition. The renormalized energy density 
is given by 

1 + AZ 



+6 



96tt 2 



-a- 4 (0') 2 + V £ (0) + 6(C + AC)a- 2 iJ0' 
(3m 4 + 12m 2 m 2 + 8aHmm' + 2a 2 # 2 m 2 ) 

r + AA +f 

<J"tt ~\ r C-fin , 



(5.29) 



K K 

where we have introduced the finite renormalizations 

AA = Km A fLf 

,2 



mi 

\Z G —K—^-Lf 



(5.30) 
(5.31) 



Though these terms have a local form they need not be considered as renor- 
malizations of Z G and of the cosmological constant, they are just finite parts 
of the mode integral. 

The divergent part of the trace T = is obtained from the one of jF div 
by multiplying with m(r)a~ A+e as 



-*fl,div — O 



,2^2 



4m 4 L + 



4tt 2 



+ 2mm" L 



(5.32) 



16 



Using the same procedure as above, and inserting all the divergent counter 
terms as they have already been determined we find the renormalized trace 

%cn = 4V ren (0) 

+6 U- a~ 2 (cC 1 ^ 1 -H$f (5.33) 

+6(C + ACK 3 U" + 2 -Ra~^) + + 4^ 

\ 3 / K K 

H (a 2 Rm 2 - 8a 2 H 2 m 2 + Y2mlm 2 + 4mm" 

+6m 4 — 2m 2 + IQaHmm'^j 

Thus far we have discussed the way of handling the infinite terms by renor- 
malization of the bare parameters. We also have considered carefully the 
finite terms left over after renormalization, in particular those that appear 
due to the continuation of the space-time dimension via divergent factors 
l/(n — 4) multiplying terms that vanish as n — 4 as n — > 4. We have con- 
sidered, however, only a restricted class of metrics, conformally flat ones. 
Furthermore, while our regularization is relativistically covariant it is not so 
under general coordinate transformations. As a consequence, the anomalous 
contributions do not appear (they neither do, for the same reasons, in Ref. 
[ PT[ ] where a Pauli-Villars regularization is used). Indeed for the Dirac field 
we seemingly need no counter terms of the form ^H^, or H^ v . In 

general, however, one expects the divergent part of T^ u for non-interacting 



fermions to be (see, e.g., |2£| |29], p0| ) 

1 r7 5 il 

L 



, im * q - 2m2 G -±\ 7 -H - h ~UH 
1 n(n — 2) 3{n — 2) 180 



T ^ [n(n-2) 9tM/ 3(n-2) G " ,,u , 

^ (5.34) 

where L has been defined in the previous section, Eq. ( |5.21| ). The first 
two terms in the parenthesis are cancelled by the cosmological constant and 
gravitational wave function renormalization counter terms, Eq. ( |5.27|) , they 
agree with the divergences we find in our analysis. The remaining terms are 
higher curvature terms, for whom we have not found any divergences. This is 
due to the fact that for conformally flat metrics in n = 4 one has the identity 

H^= (2) H^ = ^H^ (5.35) 
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and the expression in square brackets, i.e., the entire higher curvature contri- 
bution vanishes, in agreement with our analysis of divergent terms. However, 
we have to supply counter terms for those contributions as well, as for metrics 
that are not conformally flat such counter terms would be required. That 
means we have to subtract these terms. Even in n = 4 they yield a finite 
contribution, the anomaly, if one continues the expressions for those tensors 
from n ^ 4 to n = 4 and takes into account the factor 2/e = 2/(n — 4) in 
the divergent factor L. The explicit expressions for these tensors for n ^ 4 
are given in Appendix A. The finite contributions to which we call 
have been evaluated on the basis of these formulae as well as ( |5.34| ) , using 
maple. We find 

T,^ an = —^—{-3HR- —H 2 R + 22H 4 -r\ . (5.37) 
" 960tt 2 I 3 J 

The anomalous part of the energy-momentum tensor is conserved separately 
from the nonanomalous part. 



6 The initial singularity 



The fluctuation integral, the energy density, and the trace of the energy- 
momentum tensor contain contributions f S i ng , £ S in g and % ing that become 
singular as r — > 0. We have discussed previously |1| how to get rid of such 
singularities by a Bogoliubov transformation. There we had assumed that 
m(0) = 0. Here m'(0) is necessarily different from zero, as a'(0) = H (0)a(0) 2 
is different from zero due to the first Friedmann equation. So the singular 
terms are different and the singularities are stronger. For the Minkowski 
geometry this has been discussed in |5JJ. From Eq. ( |5.10| ) we have, using 
explicit expressions given in Appendix B of Ref . |34[ , 



•^sing (t) 



m'(0) , m"(0) , 



4vr 2 



471-2 



The energy density stays finite as r — > and the trace behaves as 



%ing( 7 ") 



'm'(O) _,, rn"(0), ) 



(6.1) 



(6.2) 
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As it is simply proportional to the fluctuation integral, it is sufficient to 
render the fluctuation integral finite. To do so one performs a Bogoliubov 
transformation 



cos((3 P)S )b PyS + sin(^ P)S ) exp i5 ^ s dl 
ps = -sin(/3 PiS )e-^4 PiS + cos(/3 P)S )Jl 



p,s 



(6.3) 
(6.4) 



and defines the modified initial state as being annihilated by the new oper- 
ators b PiS and dL ps . The fluctuation integral then takes the form 



E 



{^p.s^p,* sin 2 /3 P)S + V_ P)S V P)S cos 2 /3 P)S 



(2tt)32^ 
+V-^ s U^ s e iSp ' a cos/3 PiS sin/? PjS 
+U p , s V. PtS e- i5p - 3 cos/3 PiS sin/3 PiS } 



(6.5) 



or, explicitly 
f(r) ■ 



m r 



m"(r) 



m 3 (r) 



+ 



m(r)m. 



o 



m ' (0) sin(2£ r) + cos(2£ r) + i^(p, r 



2(^o) 2 
+ sin(2/3 p . 



4(S ) 



Eo 4(£ ) 3 2(^)3 2(£ )3 

(6.6) 



ps 



-21m (/ p /;e ! ^)+2mRe(/;e tf - 



2£ (£o + mo) 

Here we have used Eq. ( |5.7| ) and U P , S U P!S = —V PtS V p ^ s as well as 

2ps 



^-p,s^p,s (^p,s^-p,s 



E + m 



(V p /; + m(r)/ 2 ) . (6.7) 



In determining the parameters /3 PtS and <5 P)S the mode functions in the last 
bracket of Eq ( |6.6|) can be replaced by their leading behavior 



21m (f p f p 
2m(r)Re (/ p 2 



v 

2 AS 



2E cos(2£ r - <W) 
2m(r) cos(2E r - S PtS ) 



(6.8) 
(6.9) 



As we discuss the limiting behavior as r — > we can replace m(r) by mo- 
We then obtain the following conditions for the cancellation of the singular 
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contributions 



m"(0) 

- cos(2^ P)S )— ^ + 2pssin2/3 p>s cos(5 p , s ) = (6.10) 



2E { 
m'(0) 



cos(2/3 PiS )— ±± + 2p S sin2/? PiS sin(5 PiS ) = (6.11) 



so that the parameters are obtained as [| 

tan5 p , s = -2^o^ (6.12) 
tan2/? p , s = -^^-Jl + tanH. (6.13) 

Here we have retained only the leading asymptotic behavior. The resulting 
modifications of the various mode function integrals can be read off from 
those of Ref. One should note that m"(0) does depend, via the equation 
of motion for 0(r), on the modified fluctuation integral .F(O) which now does 
not vanish as r — > 0, and this again depends on m"(0). So m"(0) has to be 
determined self-consistently. 

The situation before removing the initial singularity is actually even worse 
than it is apparent in the analysis of ^ing. By the second Friedmann equa- 
tion a singularity of T(r) entails a singular behaviour of a"(r), but m(r) 
contains a contribution a"(r)m/, so m"(0) does not even exist. This just 
reflects the impossibility of starting the evolution consistently without get- 
ting rid, beforehand, of the initial singularity. An alternative way of avoiding 
these initial singularities has been discussed recently for scalar fields, using 



adiabatic subtraction fll"5|. 



7 Conclusions 

We have derived here the renormalized equations of motion for a massive 
Dirac field, with scalar coupling to a Yukawa field, in a flat FRW universe. 
We have used a formalism - essentially a resolvent expansion - that is useful 
for numerical computations and allows at the same time for a straightforward 

4 The right hand side of Eq. (V.9) in should have a factor 4 in the denominator, 
instead of 8, in accordance with Eq. (V.10). 



20 



analysis of the divergent terms and their separation from the finite parts 
to be computed numerically. Straightforward here does not mean simple. 
For a massive Dirac field the Feynman diagrams for the effective action are 
divergent up to diagrams with 4 external scalar fields. This reflects itself in 
a rather lengthy analysis, in our formalism as well. We have not reproduced 
all the details here, but rely on our previous work [18| for a Yukawa theory 
in Minkowski space. 

The divergent terms obtained in our formalism have been found in con- 
sistency with standard results, and the renormalization counter terms can be 
chosen independent of the initial conditions. The fact that we have consid- 
ered a massive fermion field leads to a larger number of couplings of the scalar 
field, necessary for renormalization. As these may be obtained by a shift of 
the scalar field, this modification is a straightforward generalization. The 
renormalization counter terms have been determined in dimensional regular- 
ization and with a MS convention for the renormalization. The formalism 
can be adapted to other regularization and renormalization schemes. 

The trace anomaly does not appear explicitly in our formalism. This is 
a consequence of the special geometry we consider here and does not imply 
an inconsistency of the formalism. This can be understood from the general 
analysis of the divergent parts of the energy-momentum tensor, as given, e.g., 
in p0| . We have used this analysis in order to infer the anomalous part of 



the energy-momentum tensor. 

We have also analyzed the initial singularities. They obstruct the cos- 
mological initial value problem in an essential way, a singularity in the trace 
of the energy momentum tensor leads to a singularity in the expansion rate. 
So they have to be avoided by a suitable choice of the initial state. We have 
constructed explicitly such a state; it is not unique but can be modified if 
only the asymptotic behavior as p —>■ oo of the Bogoliubov parameters is 
retained. 

Our results can be used in a straightforward way for numerical compu- 
tations, which we plan to carry out in the near future. We do not expect 
an essential modification of results related to the parametric resonance phe- 
nomenon; however, especially for large Yukawa couplings and/or masses, 
the finite remnants of the renormalization counter terms and the singular 
parts are not necessarily small and this signals the importance of considering 
renormalization carefully. 
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A The higher curvature tensors 

The higher curvature tensors are defined by 
1 5 



l <) 



Jd n x^R 2 



iw ~~ %9 iiv °R — -^9 + 2RR n V , 



(A.l) 



2i*« -OR 



Jd n x^gR aP R 



ai3 



H 



1 5 



,"" " ~j9 W D ^ + ^ au ~ 1^9 nvR R a/3 



fJbV 



Jd n x^R a ^ s R 



af3-y8 



- l -g^R a ^ s R aPl& + 2R"^"Rf^ - AUR^ + 2R 
AR m RZ + 4R a(3 R 



[IV 



afj,f3u 



(A.2) 



On order to evaluate the anomalous contribution to the energy-momentum 
tensor we need the continuation of the higher order curvature tensors for the 
FRW geometry to n = 4 — e dimensions. We choose to continue the spatial 
dimension to 3 — e. We already have given the n- dimensional continuation 
of the Ricci tensor in Eq. (|2.7| ) . 

For the time-time components and the trace of the tensors ^ n >H^ v we 
obtain 



W H tt 



-6HR + ^R 2 - 6H 2 R 



(A.3) 
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+ {n - 4) (-2HR - {n + l)RH 2 ) , 
{2) H tt = -2HR + ^R 2 - 2H 2 R + (n - 4) {^\hR (A.4) 

-24(^1) -j(» + 2)^ + i( B -l)( B -2)^) , 
# tt = -2HR+^R 2 -2H 2 R (A.5) 

£ = -6i2 - 18HR + (n - 4) (-2i2 - 2(n + 2)/LR - ^i? 2 ) (A.6) 
( 2 )#£ = -2R-6HR+ (n-4) - ^(n + 3)iLR (A.7) 



2 2 
ni? 2 1, 1 



+ ^(n - 2) 2 # 2 i? - -n(n - \){n - 2) 2 H* 



8(ra-l) 4 V y 8 

(p2 
-2M - - 
z(n — ij 

+ {n-2)H 2 R-^n{n-l){n-2)H A ^ . (A.8) 

B The Dirac equation in n dimensional FRW 
geometry 

The Dirac equation in curved space is formulated using the n-bein e a = e^dx 11 
and the spin connection dx u . It has the form 

i>fK [d v + cof^ [ 7a , 7*]) i> - rmf; = . (B.l) 

Here Eg is the inverse n-Bein and the connection is defined by the equation 

dAe a + u%Ae b = (B.2) 
Choosing for the FRW geometry e° = dt and e k = a dx k one finds 

u; k o = ddx k , (B.3) 
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while the purely spatial components u k j vanish. Using E® = 1, E\ = 1/a, 
all other components being zero, we obtain 

i (l% + - 7 fe V fe + - 7 ^ feo i [ 7fe , 70]) ^ - mi, = , (B.4) 
\ a a 4 / 

or, using 7 fc [ 7fe , 7o ] = 2(ra - 1) 7 ° 

^7° (do + ^ + ( Vv fc - m) ^ = . (B.5) 

It is the almost obvious generalization of the Dirac equation in the n = 4 
FRW geometry. 
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